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We study the build up of quantum coherence between two 
Bose-Einstein condensates which are initially in mixed states. 
We consider in detail the two cases where each condensate is 
initially in a thermal or a Poisson distribution of atom num- 
ber. Although initially there is no relative phase between the 
condensates, a sequence of spatial atom detections produces 
an interference pattern with arbitrary but fixed relative phase. 
The visibility of this interference pattern is close to one for the 
Poisson distribution of two condensates with equal counting 
rates but it becomes a stochastic variable in the thermal case, 
where the visibility will vary from run to run around an aver- 
age visibility of 7r/4. In both cases, the variance of the phase 
distribution is inversely proportional to the number of atom 
detections in the regime where this number is large compared 
to one but small compared with the total number of atoms in 
the condensates. 

PACS numbers: 03.75 Fi, 05.30 -d 



I. INTRODUCTION 

The recent experimental realization of weakly inter- 
acting Bose-Einstein condensates has stimulated a 
large amount of theoretical work on the properties of 
these condensates. Recently, there has been a great inter- 
est in the interference and the establishment of a relative 
phase P-pd[ between Bose-Einstein condensates which 
start in Fock states or as a mixture of coherent states. 
In this paper we show by numerical and analytical calcu- 
lations how a relative phase is established between two 
independent condensates initially in mixed states. In par- 
ticular we look at the quantum interference when the 
condensates begin in thermal or Poisson mixed states. 
There is a significant difference in the visibility of the 
interference patterns between these two cases. 



II. QUANTUM PHASE BETWEEN INITIAL 
MIXTURES 

We consider two Bose-Einstein condensates which are 
dropped on top of each other. This system was first pro- 
posed by Javanainen and Yoo M. An example of such a 
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system is a hot wire grid placed below the condensates 
which acts as an atom detector since it removes atoms 
from the condensates as they fall under the influence of 
gravity. Javanainen and Yoo showed that a spatial inter- 
ference pattern would be observed illustrating the pres- 
ence of a relative phase between the falling condensates. 
We will show here that this relative phase is also present 
when the initial states of the condensates are no longer 
a pure Fock state of known number but are a mixture 
of number states. This corresponds physically to the sit- 
uation where we are uncertain of the initial number of 
atoms in each condensate. We shall take this mixture to 
have a number distribution P n giving an initial density 
operator of the form 
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oo 



P= ^ 

n\— ri2— 



P ni \ni)(ni\ <g> P n2 \n 2 )(ri2\ 
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where n\ and n 2 refer to the first and second condensate 
respectively. The spatial interference is established since 
when we detect the atoms from the condensates we are 
unaware which condensate the atom came from. This 
interference pattern is built up from the spatial detec- 
tions of individual atoms. After m detections we observe 
atoms at positions {x\, . . . , x m }. If we denote the field 
operator for the detection of an atom at x as ip (x) then 
the joint probability of m detections at the {x\, . . . , x m } 
positions is 



P 



(m 



\ Xl , ...,x m )= AA (m) Tr [p^{ Xl ) . . . ft(x m ) 

x$(x m )...'ip(xi) \ , (2) 



where the symbol Tr denotes the trace over the n\ and 
ni number states. The normalization A/"' m ' is defined by 



dx\ . . . dx m Tr 



•Em ) 



x jp{x ri 



(3) 



which is independent of x\ , . . . , x m but will in general 
depend on p. We define the field operator for the two 
condensate system pLfLOl as 



ip(x) = &i + VF 



a 2 e 



-i(f>(x) 



with 



= (ki — k 2 )x 



(4) 
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where a\ and a 2 are the atom annihilation operators for 
the first and second condensate respectively with each 
condensate possessing momentum k\ and k 2 . We define 
the ratio T — 72/71 where 71 and 72 are the detection 
rates for each of the condensates. We shall assume 71 > 
72 without restriction of generality. 

Substituting equation (|l|) into (||), we obtain the fol- 
lowing expression for the joint probability. 

OO OO 

= A f(m) J- J2 Pn u n 2 {ni,n 2 \^{x 1 )..4H.Xm) 
rii— n 2 — 

xip(x m ) . . .^(xi)\n ll n 2 ) (6) 

where we have written P ni ,n 2 = PnxPn^ for convenience. 
The expression for the joint probability above shows that 
it is a weighted sum of probabilities over fixed numbers 
n\ and n 2 of initial numbers of atoms in each conden- 
sate. The weighting is determined by the probability 
distribution of the number of atoms. For very narrow 
distributions, corresponding to an initial fixed number of 
atoms in each condensate, we see that this will simplify 
down to the expression used previously by Javanainen 
and Yoo [Q. However for broader distributions, for ex- 
ample a thermal distribution, this sum will affect this 
probability and hence the spatial interference. 



A. Visibility conditioned on 1 detection 

Before we show how to numerically generate this inter- 
ference pattern, we look at the build up of interference 
for the first two detections. Let us start by considering 
a joint Fock state. After one detection in which an atom 
is observed at x\, the un-normalized state vector for the 
system is 



+ Jf^' 2 e- l ^ Xl) 



\ni,n 2 - 1). 



(J) 



The joint probability density for detecting atoms at x 
and X\ starting from \ri\, n 2 ) is then proportional to 

(cp 2 (4>) \<P2 (4>)) = n 2 |^ t (xi)-0 t (a;)^(a:)-i/'(xi)|ni, n 2 ) 
= nx(nx - 1) + r 2 n 2 (ri 2 - 1) 

+2T ni n 2 {l + cos[(j){x)- 0(xi)]} . (8) 

Since the initial state of interest is a mixture of Fock 
states with weights P nii „ 2 , the joint probability for this 
initial state by Eq. (H) is 



P 



(2) 



(x, Xl ) = AA 2 ) 



P„ 



(n 1) n 2 \'ii)^(x 1 )^(x) 



xip(x)'il>{xi)\n 1 ,n 2 } 



2) 



Pn 



(ni(ni - 1) + T 2 n 2 (n 2 - 1) 



+2r ni n 2 {1 + cos [4>{x) - 0(xi)]}) 
= ^ 2 )([{n?)-(n 1 )]+r 2 [(nl)-(n 2 )] 

+2T( ni )(n 2 ) {1 + cos [(j){x) - {x x )]}) (9) 

where we have used angle brackets to denote averages 
taken over P ni i7J , 2 . These factorize as we assume that 
Pni,n 2 — PnxPn?.- The conditional probability y(x\x\) — 
{x,x\) /p^ {x\) differs from the above only by a x- 
independent factor and we may write 



P(x\xi) = ^^-{1 + Vcos [4>(x) - <j>(xi)]} 



where 
V = 



2Y{ ni ){n 2 ) 



[<n?)-<m)]+ra [(ni)-(n 2 )] + 2r(n 1 )(n 2 ) 



(10) 



(11) 



may be interpreted as the conditional visibility of the 
interference pattern. 

For a thermal distribution we use the following re- 
lationship between the second moment and the mean 

TO = (to), 



(to 2 ) = 2to 2 



(12) 



Substituting the above equations into Eq. ([11]) we obtain 
an expression for the conditional visibility in terms of the 
means 



Vthcrmal (ni,n 2 ) ■ 



Tni n 2 



n\ + (Tn 2 y + Tn 1 n 2 



(13) 



which gives a maximum visibility for equal net detection 
rate in each condensate, i.e. giving a maximum value of 
one third for t?i = Tr^. 

Alternatively, for the case of a Poissonian number dis- 
tribution we use the relationship 



(m 2 ) = TO 2 + TO. 



(14) 



Proceeding in the same manner as in the thermal case 
we obtain 



Vp 



i(ni,n 2 ) 



2Tn\ ri2 
(ni + Tn 2 f 



(15) 



which has a maximum value of one half for TL\ = r?^. 
The values of one third and one half that we have ob- 
tained have also been seen in optical experiments where 
intensity correlations were measured for Poissonian and 
thermal light sources by Rarity et. al. jjjj 

Let us also consider the limiting case of fixed initial 
number, where we know that there are exactly n\ atoms 
in one condensate and n 2 in the other. Then using Eq. 
( |i"l| ) the conditional visibility is 



VFock(«l,n2) 



2rni n 2 



(rii + Tn 2 ) 2 — (ni + T 2 n 2 ) 



(16) 



2 



The maximum occurs again when n\ = Tn 2 where the 
net detection rates for both condensates are equal. This 
is not surprising since the size of the interference term 
depends on our lack of knowledge of which condensate a 
given detected atom comes from. If the net detection rate 
of one condensate is larger than the other, then we know 
that the detected atom is more likely to have originated 
from this particular condensate. Unlike the previous two 
cases, the maximum visibility depends on TV the initial 
total number of atoms in the condensates. The maxi- 
mum visibility for the special case of equal number and 
detection rate for each condensate we have 



[V, 
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Fock. 
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2(1 - l/N) 



(17) 



In the limit of large N the conditional visibility ap- 
proaches that of the Poissonian mixture which has the 
value of one half. For practical purposes where the num- 
ber of atoms in each condensate is well over one thou- 
sand, the cases of a Poissonian mixture and of initial 
Fock states are indistinguishable. We plot these condi- 
tional visibilities for the different initial conditions as a 
function of the ratio of initial net detection rates, Tfii jn 2 
between the two condensates in Fig. ([!]). To see how the 
maximum visibility changes as we alter the width of these 
initial distributions, let us also consider some arbitrary 
Gaussian distribution with a variance a 2 and mean n. 
The conditional visibility when both condensates start 
with this Gaussian distribution and equal detection rates 
is 



V G 



aussian 



(n,a 2 ) 



2n z 



(18) 



This is approximately equal to (2 + cr 2 /n 2 ) 1 for large 
n. For wide distributions a 2 ^S> n 2 , we see that the con- 
ditional visibility tends to zero. Conversely for narrow 
distributions where a 2 <C n 2 , the visibility becomes ap- 
proximately one half for large n. In the special case of 
(7 2 = n where we approximate the Poisson distribution by 
a Gaussian, Eq. ( [l8| ) yields a visibility of one half which 
is consistent with the value obtained from the expression 
for the Poissonian visibility given by Eq. (|lfj) . 



B. Conditional probability density after m detections 

The conditional probability density p(x\xi) displayed 
by Eq. ( fL0| ) can be generalized to an expression govern- 
ing the probability density of x given the previous m 
measurements {x\...x m }. We can write the operator 
representing the cumulative effect of m detections as 



* (x m ) * (x m _i) . . . <l (xi) = Y[ (ai + VTa 2 t 



3 = 1 



(19) 
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where we define 4>k = <p{xk) for notational convenience 



and the coefficients 7r 



(m) 



j) can be found by 



computing the power series expansion 
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,4> m )z k . (21) 



They satisfy the recursion relation 

^+1) = ^(ro) (1 _ 6k m+i) + ^) (1 _ Sk o) g - #m+1 

(22) 



(rn) _ (m) 



») for 



where 

we have used the notation 7r 
brevity. In a numerical simulation, the product (|2l| ) can 
be updated after every atomic detection by carrying out 
polynomial multiplication. The un-normalized state vec- 
tor after applying the above operator to the initial state 
\ni,n 2 ) is 



|<Pm) = 



ni\n 2 \ 



(ni — - m + fc)! (n 2 — fc)! 



X7r (m) r fe/ 2 | ni _ m + k,n 2 - k). 



(23) 



Let us now consider the un-normalized wave-function af- 
ter the (m + l)'th detection 



\(p m +i (<f>)) = (&1 + Vfa 2 e \ip m ) 



(24) 



where we have explicitly shown the <fi dependence. The 
joint probability of m+1 detections at the {xi, . . . , x mi x} 
positions is 

^ (^1 ) • • • 3 X m , x) 



: AA(" 1+1 )Tr [^{ Xl ) . . . ^(x m )^(x) 
x^(x)ip(x m )...tp(xi)^ (25) 
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:Af(m+ l) J2 J2 PnunAfm+1 (<t>) Wm+l (</»)), (26) 
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where 



(tpm+i (</>) \<Pm+i ((/>)) 

m — 1 
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r fe+1 ni!n 2 ! 



(ni — m — 1)! 



k=0 



(ni — m + fc)! (n 2 — k — 1)! 



r m+1 n 2 ! 



(n 2 ~ m - 1)! 
The conditional probability density is thus 



(27) 



fe=0 
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p(x\xi, . . .,x m ) 

- TVM [{ ni ( ni _ l) . . . ( ni _ m )) 



+ 



m— 1 

E 

fc=0 



(nx (ni — 1) . . . (ni — m + k + 1)) 



x(n 2 (n 2 -l)...(n 2 -fc))r fc+1 
+r m+1 (n 2 (n 2 - 1) . . . (n 2 - m))] 
where 



p(™) (a?i,..., z m ) 



(28) 



(29) 



is a a; independent normalization factor. The angle brack- 
ets denotes the sum over the probability distribution 
p 

1 III ,™2 • 

For different initial mixtures, either a thermal or a 
Poissonian distribution, we obtain different relationships 
between the higher order moments and the first (i.e. the 
mean). The relationship for the thermal case is 



(n (n - 1) . . . (n - fe))thermai = (* + 1 )! n k+1 (30) 



and 



(n (n - 1) . . . (n - k)) Po 



-k+l 



(31) 



for the Poissonian distribution. Using these properties we 
obtain the conditional probability density for the thermal 
distribution 



PPo 



(x\ Xl ,...,x m ) =N^l son [nT +1 + (Tn 2 r +1 
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(Poisson) 
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(32) 



and 



Pthcrmal {x\x X , . . . , X m ) = ^t^mal { ( m + X ) ! K™ +1 

m. — 1 

+ (rn 2 ) m+1 j + ]T .4 



(thermal) 



fe=0 



(m) 
fe+1 



(33) 



for the Poissonian distribution. Note that the x inde- 
pendent normalization factors are now calculated with 
respect to the relevant density operator, i.e. Poisson 

for ^PoiLon and thermal for N^ Imal , Both the thermal 
and Poisson distribution's conditional probability densi- 
ties have a similar form with different weighting factors 
A k given by 



A 



(Poisson) 



- rn — k — k + 1 k + 1 

n-. n J- 



(34) 



and 



A 



(thermal) 



(to - fc)! {k + iy.A 



(Po 



(35) 



The generalized conditional probability distributions Eq. 
( |32| ) and Eq. (|3^) groups all the measurement depen- 
dent terms, the {4>i, ■ ■ . , 4> m } detections, within the |. . .| 
brackets. Numerical calculations are readily obtainable 
since we need only to simulate the expression within the 
|...| brackets. We shall come back to this point and 
explain in more detail the numerical simulation proce ss 
when we describe the numerical results in section [I D . 



C. Analytical results for a large number of 
detections 

1. Poissonian mixtures 

In the case of Poissonian mixtures the expression (|3J 
can be made more explicit and analytical results can be 
extracted. As shown by Cirac et. al. [(7) it is useful to 
represent the state after m detections in the form of a 
P-representation 



with 



n\, a 2 



'n 2 . 



(36) 



(37) 



In writing Eq. ( p6[ ) we are assuming that the observation 
time t m satisfies "fit m <C 1, so that only a negligible frac- 
tion of the atoms in the condensates are counted. The 
initial Poissonian mixture po is also of the form of Eq. 
© with 



/o W = I- 



(38) 



The effect of an additional counting event on ( pq ) is to 
change 



with 

fm+i (VO 
where 



N, 



fm -> fm+1 (tp) 



[l + Acos(0 m+ i-V)]/m(V') (39) 



m+1 



2y/Tnin 2 
fix + Tn 2 



(40) 



and N m+ i is determined by normalizing J -^f m +i (VO = 
I after each counting event. We note that < A < I. 
In the following we shall assume that A < I, as the case 
A = f requires a separate mathematical treatment. The 
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physical meaning of f m (ip) is clear from Eq. (|36|): it is 
the probability distribution of the relative phase between 
the two condensate modes. The explicit form of fm (ip) 
as a function of ip is easily obtained from Eq. ( |38|) and 
Eq. © asQ 



m 1 

f m W = Y[[l + \cos(cf )k ~^)}—. 



k=l 



(41) 



The normalized joint probability distribution to observe 
the phases fa,..., <p m +i is equal to 



/, , m+l m+l 
fc=i fc=i 



(42) 



Therefore the conditional probability ( 32 ) determined by 
the ratio of 2 joint probabilities is equal to N m +i and 
given by 

/dip 
— [1 + A COS (<p - Ip)] f m (ip) ■ 
Zir 

(43) 

where </> = </> (x) and f m (ip) depends on the previous m 
detections. Now we investigate the behavior of f m (ip) 
as a function of ip for large m. We shall assume that 
it becomes a narrow distribution centered around some 
m-dependent maximum ip m with a variance a m , i.e. we 
put 



fm(lp) 



1 



/ 27TCT T 



■ exp 



(ip - 1p m Y 

2al 



(44) 



We assume here that a m < 1, so that the 27r-periodicity 
of fm (ip) is not in noticeable conflict with (f44|). 

We shall now determine the time evolution of ip m and 
af n , using Eq. (|39|), and show that our assumptions are 
self-consistent in that indeed a m <S 1 for m ^> 1. From 
Eq. ( |4l| ) we obtain by taking the logarithm 



In fm (ip) 



fe=i 



In [1 + A cos (<pk — ip)]+ c 



(45) 



where the constant c depends on the (p m but not on ip. 
The maximum ip m of In / m (i/>) must therefore satisfy 



A sin (fa - ip m ) 

- 1pm) 



' 1 + A cos ( 

k— 1 



= 0. 



(46) 



x This formula permits us to appreciate wh y t he case A = 1 
is very special: For A < 1 the distribution (Wll) is a positive 
function f m (ip) > everywhere. For A = 1 it has m 2- 
fold degenerate zeros, i.e. f m (ip) cannot approach a smooth 
function for m — > oo. 



The evolution of the maximum is obtained similarly by 
taking the logarithm of Eq. ( ^) and using Eq. (|44|) . We 
find 



(i'm+l ~ 1pm) 



A sin ( 



"m+l 



1pm) 



'm+l 



1 + A COS ( 



"m+l 



1pm) 



= 0. (47) 



ft is clear that (ip m +\ — ipm) depends on the outcome 
of the (m + l)'th measurement. If f m (ip) is, in fact, a 
narrow distribution on the scale 27r, as we have assumed, 
then the probability distribution to find <p>m+i in that 
measurement can be estimated as 



P, 



m+l (<Pm+l 



) ~ 1 + Ae~ 



l+ i - ipm) (48) 



as follows from Eq. ( ]43| ) and Eq. ([H|). Based on this we 
find 



(ip 



1pm+l - Iprr 
2 

'm+l 

IpmY 



= o, 



m+l 



a 



m+l 



2n 1 + A cos ( 



d<p cos <p sin <f> 
2tt(1 + Acos0) : 



= 1- Vl -A 2 

1 2 (1- VI - A 2 ") 2 , 

vt-a^ (49) 

where we used <r m C 1 to expand to first order 
exp (—am/2) = 1 — of rl /2. We conclude that, if A is 
bounded away from 1 , for small variance of the phase dis- 
tribution the variance in the jitter of the average phase 
ip m in subsequent measurements is of the order of the 
square of the variance a m in the phase distribution, i.e. 
the position of the maximum is very stable and may be 
considered as fixed in the limit we consider. 

Turning to the time evolution of the variance we ob- 
tain again using Eq. (Q) on the right hand side of Eq. 
( p9| ) , taking the logarithm and the second derivative with 
respect to ip 



1 



'm+l 



A 2 + A COS ((p. 



m+l 



Ipr, 



o m [1 + A cos ( 



- Ipn 



(50) 



ft shows that typically the inverse variance grows accord- 
ing to 



1 



1 



'm+l 



0(1) 



(51) 



which is the reason why the variance itself indeed be- 
comes small. Averaging as before we obtain for A 
bounded away from 1 



1 



1 — \/l — A ; 



O(al) 



const. (52) 
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which shows more explicitly how the inverse variance 
grows on the average and becomes large for large m. 

Finally, we relate the phase distribution f m (ip) after 
m measurements to the observed interference pattern. 
Experimentally, a phase distribution may be extracted 
by fitting the observed interference pattern, normalized 
with respect to its constant part, to the expected density 
Eq. (pl3| ) after m counts 

/^[l + Acos^-VOl/mWO = l + A'cos(0-^ m ) 

(53) 



with 



A' = Ae" 



(54) 



Here we used the Gaussian form of f m (ip) assumed in Eq. 
(fl4|). Equations ((5^), (|54]) allow us to extract numbers 
for ip m and <r^. A weak link in the argument leading 
to Eq. ( f53| ) might seem to be the fact that the condi- 
tional probability (43) is proportional to the expectation 
value of the density in the condensates, while what we 
really need is a normalized measure of the density of the 
counted atoms. However, all that is really required is 
that these two quantities should be proportional to each 
other, which is an assumption implicitly already made 
when assuming that the counting rate is proportional to 
the number operator in the condensate. 



2. Thermal mixtures 

If the initial states of the two condensates are thermal 
mixtures we may adapt the results for the Poissonian 
mixtures as follows. The thermal initial state can be 
represented as a mixture of Poissonian states via 
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Po 



Z1Z2 

ni ,n 2 



V- e -/J(n 1+ na) L> /„. 



= -j-j. I d 2 a\ I d a% exp 



fix 



n 2 n 2 
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\ai) 



\ a 2 I 



Jo 



dxi exp 



X2_ 



[dip f dip! 



(55) 



with the inverse temperature proportional to f3 and nor- 
malization factors AT and M' . The evolution of this initial 
state due to the counting of atoms can now be obtained 
by using the evolution for Poissonian mixtures for fixed 
Poissonian averages x\ , X2 under the integrals over x\,x%. 
We obtain with / (ip) — 1, using Eq. ( j4l| ) for 



A = A (xi,x 2 ) 



2y/rxxX2 
xi + Tx2 



(56) 



771 

x 2 cos (4> k — ip) \ ■ (57) 



k=l 



The optimal visibility A thus varies according to Eq. (|56| ) 
for different members of the ensemble whose x\, x 2 values 
are exponentially distributed according to Eq. (pq) . The 
average visibility becomes 



A = 



(xi)+r(a;2> 



= ^ 2y/r ( Xl ) (X 2 ) 
4 (x 1 )+T(x 2 ) 

which may be written asn 



(58) 



For equal average counting rates 
fix = Tn 2 

the average visibility for initially thermal mixtures is sim- 
ply 



A 



4' 



(59) 



which agrees well with the result of the numerical simu- 
lations to be presented below. 



D. Numerical results 

The form of the expressions for the cond ition al prob- 
abilities, Eq. (32) and Eq. ( |33| ) of section (ilB) can be 
readily applied to stochastic simulations of atom detec- 
tion which generates interference patterns. The proce- 
dure follows the spirit of Javanainen and Yoo's work; a 
random number (f>i is generated for the initial atom de- 
tection (since we know the initial conditional probability 
distribution is uniform), then it is used to calculate the 
conditional probability density p{<p\4>\) so that the sec- 
ond detection <p2 can be generated by randomly selecting 
a value according to p(<p\(f)i). This process is repeated 
to generate m atom detections which are binned and dis- 
played as an interference pattern. 

Examples of these stochastically generated interference 
patterns are displayed in fig. ^| as histograms of the raw 
output from the simulations, the sequence of detected 



2 This is consistent with the visibility defined by Eq. ( flc| ) 
previously since A ((xi) , (^2)) = A. 
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positions {4>i, ■ ■ 4> m } are sorted into 25 bins plotted as 
circles. Figure |2|(a) shows the interference pattern for 
an initial thermal mixture which has a visibility of 0.79 
whereas in fig.|||(b) we have a visibility of 0.97 for the 
Poissonian mixture. These visibilities are calculated via 
a least-squares fit of the form 1 + (3 cos (2tt + </)) shown 
as the solid curves in fig. ||(a) and (b). In both cases we 
have simulated 500 detections. 

The numerical simulations used to generate the his- 
tograms in fig.^(a) and (b) also calculate the conditional 
probability distributions before each detection as these 
are necessary to calculate the location of the detected 
atom. The evolution of these conditional distributions 
gives an insight into the build up of the interference 
pattern juj. The visibility of the complete interference 
pattern (as calculated from the least-square fit) can be 
considered as an average over these conditional visibil- 
ities. Graphs of the conditional visibility as a function 
of atom detections typically approach a value close to 1 
within 100 detections and stay at that value thereafter. 
Thus the conditional visibility after m detections can be 
thought of as a good approximation of the visibility of 
the complete interference pattern for values of m > 100. 
We have not displayed these graphs of stochastically gen- 
erated sequences since they possess fluctuations about a 
generic shape which is a property of all such sequences. 
We plot the average conditional visibility over many such 
sequences in fig. |] and fig. [| for Poisson and thermal 
mixtures respectively. The ratio T^/fii was set at 1 with 
the average performed over 1000 runs (we shall hereby 
refer to an individual sequence of detections as a single 
"run" for convenience). The shaded regions around the 
average conditional visibility displayed as a solid line de- 
pict the extent of the fluctuations for an individual run. 
The boundaries of this shaded region corresponds to the 
upper and lower quartiles respectively, 25% of the data 
lies below the lower quartile whereas 25% lies above the 
upper quartile, thus the probability that the fluctuations 
lie within this region is 50%. These fluctuations are much 
larger for the thermal mixtures in comparison with those 
of the Poissonian mixtures. 

The maximum value of the visibility in fig. || occurs 
at 500 atoms detected with a value of 0.999 (3 sig. fig.). 
The value of the visibility obtained from Eq. (|54|) using 
°m ~ V m i s 0.999 (3 sig. fig.) which agrees very well 
with the simulation. In the case of the thermal mixture, 
fig. |[ the maximum value of the visibility is 0.777 again 
at 500 detections. This is close to the analytical value of 
tt/4 fa 0.785 predicted by Eq. @. 

The large difference in the fluctuations between the 
two cases is due to their differing degree of sensitivity 
to particular runs. In the Poisson case, it is relatively 
insensitive since the majority of the terms in the |...| 
brackets of Eq. ( |3^ ) depends on the {<pi, . . . ,<p m } detec- 
tions fairly equally whereas in the thermal case we see 
there is an additional factorial factor in front of the | . . . | 
brackets which favors k values at the ends (close to zero 
and m— 1). This enhances the sensitivity of the visibility 



upon particular combinations of 4>k thus the sensitivity is 
high. Alternatively, the large fluctuations of the visibility 
in the thermal case is not surpr ising when we consider the 
analytic treatment of section ( II C 2) . The thermal state 
was represented as a mixture of Poissonian states so that 
the averaged conditional visibility has an additional av- 
erage over the mixture of Poissonian states. 

The variance for the Poisson case can be estimated 



from fig. 
small. W 



by using Eq. ( p4| ) since the fluctuations are 
len the variance are small compared to 1, the 



variance is 



a 2 m = 2 (A - A') 



(60) 



where A' is the numerical averaged visibility and in the 
case of equal counting rates A = 1. The variance for the 
Poisson case is graphed as the solid curve in fig. |B| with 
the dashed curve displaying the cr^ ~ 1/m relationship 
predicted by Eq. (|52|). Note that we have not plotted 
variances below 50 detections since Eq. (^4|) is only valid 
for number of detections m ^> 1 . As expected, the agree- 
ment between the two curves becomes better as more 
atoms are detected. We cannot obtain a good estimate 
of the variance in the same manner for the thermal case 
because the fluctuations in the visibility are large. 

We will show in the appendix that the results for mea- 
surements induced phase distribution for initial number 
states are the same as those obtained for the Poissonian 
mixtures. Numerical simulations of the time evolution 
and effect of detections on the wave-function of the con- 
densates have been performed jl(| where the variance 
can be calculated directly from this wave-function. This 
is possible in the simpler case of initial number states. 
Because the phase is periodic it is convenient to use the 
following measure of the spread of the phase distribution: 



6<p = l- (cosA</>) 2 - (sinA<^) 2 

where the trigonometric operators are defined as 
1 



cosA0= - 



and 



sinAfli = — ( e 



(61) 



(62) 



(63) 



with the SG phase operators e iA ^ defined for the rela- 
tive phase between the two condensates. This measure 
ranges from zero to one with values close to zero agreeing 
well with the actual variance. Since our expressions are 
only valid for small variances, this measure is very use- 
ful as an estimate of the variance. Figure || plots three 
curves of 6<j), each one with differing relative count rates 
between the condensates (TN 2 /Ni) where N-y and N 2 are 
the initial atom number in the first and second conden- 
sate respectively. Analytical predictions of the gradient 
for each relative count rate are obtained from Eq. (|5^). 
The solid, dashed and dash-dotted curves of fig. || cor- 
responds to TN2/Ni ratios of 1, 1/2, and 1/4 which are 
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predicted to have gradients 1, 2/3, and 2/5 respectively. 
The numerical gradient was calculated for points from 
20 detections onwards, the points below 20 were ignored 
because the predictions required to 3> 1. The curves 
clearly show larger curvature in this region in compari- 
son to points after 20. We obtained numerical gradients 
of 0.97, 0.66, and 0.41 for the TN 2 /Ni ratios; 1, 1/2, and 
1/4. This good agreement between the analytic and nu- 
merical gradients verifies the relationshi p o f th e vi sibility 
A with the relative counting rates, Eq. ( |40|j5^ , |A^ ) . Note 
that if we had graphed 5<f> itself instead of its inverse, we 
would had seen a curve shaped like that of fig. hence 
the straightness of the curves of fig. ^ indicates th e acc u- 
racy of the 1 /to prediction of the previous section (|II Cffl 



III. SUMMARY 

We have analyzed in detail the build up in quantum 
coherence between two Bose-Einstein condensates which 
are initially in a thermal or Poisson state. Interference 
patterns are produced via spatial atom detections which 
establishes an arbitrary but fixed relative phase between 
the condensates. In the regime where the total num- 
ber of atoms detected is only a negligible fraction of the 
atoms in the condensates although the actual number 
of detections is much greater than 1, we find that the 
visibility of the interference pattern for the Poisson dis- 
tribution depends on the relative counting rates of each 
condensate with a maximum of one for equal rates. In the 
thermal case, the visibility becomes a stochastic variable 
which varies from run to run around an averaged value 
determined again by the relative counting rates but the 
average has a maximum of 7r/4 for equal rates in good 
agreement with our numerical simulations. The inverse 
variance of the phase distribution grows linearly with 
the number of detected atoms to in the regime where 
1 <C to <C iVi, N 2 . This has been shown analytically for 
the Poisson state and this relationship also holds true for 
the thermal state since we may write the thermal state 
as a mixture of Poisson states. In the appendix we have 
shown that the results for the initial number states fol- 
lows those derived for the Poisson state. In particular, 
the inverse variance is proportional to m, which has been 
numerically verified. 

This research was supported by the University of Auck- 
land Research Committee, the New Zealand Lottery 
Grants Board and the Marsden Fund of the Royal So- 
ciety of New Zealand. 



3 Note that section (II C ) consider Poisson and thermal states 
but fig. [] show results from simulation of initial number 
states, however we will show in the appendix that the width 
of the phase distribution has the same relationship with the 
number of detections m as the Poisson state in the limit 
1 < m < JVi.iVa. 



APPENDIX A: INITIAL NUMBER STATES 

Let us consider the case of initial number states and 
compare the results with those of initial Poissonian mix- 
tures. The initial number state 



|Vo) = |JVi) |JV a ) 



(Al) 



evolves into 

m 

l^m) = Y, C n{m)\N 1 -n)\N 2 -m + n) (A2) 
n=0 

where C n (0) = 5 n $. We shall assume 7iATj > 7 2 7V 2 and 
that the parameter 



2vrjVi7v 2 



(A3) 



Ni + TN 2 

is bounded away from 1. We shall only consider the case 

l<m«JVi,JV 2 . (A4) 

Then the evolution of the C n (to) under the atom count- 
ing process is 

C n (to + 1) ~ ^/n[ (1 - $„ „) Cn-i (to) 



+ v /nV2(l-^ m+1 )e-^+ 1 C„(TO). (A5) 

We introduce a phase-representation by the Fourier 
transform 

m 

Fm M = C n M e ~ mV 
n=0 

C n (to) = J ^e m *F m (ip) . (A6) 

The number density after to measurements is then 

P m (x) ~ (ip m \ ft (x) $ (x) \ip m ) 

{m 
51 (i - a»,o) 
n=0 



xe 



:[ip-4>( 



which for to 3> 1 can be approximated by 

P m (x) ~ 1 + Ajv J J \F m (v?)| 2 cos [<p-4> (x)] . (A7) 

This formula should be compared with Eq. ( flS) ) for Pois- 
sonian mixtures. It can then be seen that \F m (ip)\ and 
f m (4>) play identical roles. The evolution of F m (ip) un- 
der the atom counting process follows from Eq. (A5) and 
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F, 



m+1 



(ip) ~ (V^e"^ + v ^iV2e- i ^+ 1 ) fo>) (A8) 



and hence, providing a normalization factor N m +i to 
keep |F m+ i ((p)\ 2 normalized if \F m+ i (ip)\ is 

\F m+1 (ip)\ 2 ~ [1 + Ajv cos (^ m+ i - </?)] |F m (<^)| 2 — - — 

(A9) 

which, to the accuracy we have considered here, is the 
same as Eq. ( p9| ) for Poissonian mixtures. Hence, the 
results for the measurement induced phase distribution 
for number states in the limit 1 -C m <C N\ , N2 are the 
same as the results for Poissonian mixtures obtained in 
section 



IIC 



FIG. 3. Plot of the conditional visibility averaged over one 
thousand runs versus the number of atomic detections for an 
initial Poisson state where ni = YH2 . The shaded region cor- 
responds to the interquartile range of the individual runs with 
the mean over all runs depicted by the solid line. Therefore 
50% of the run lies within this shaded region. 



FIG. 4. Plot of the conditional visibility averaged over one 
thousand runs versus the number of atomic detections for an 
initial thermal state where ni — Yn,2. The shaded region cor- 
responds to the interquartile range of the individual runs with 
the mean over all runs depicted by the solid line. Therefore 
50% of the run lies within this shaded region. We also plot a 
dashed line at 7r/4 corresponding to the value of the average 
visibility predicted from the analytical work. 



FIG. 5. Plot of variance versus the number of detected 
atoms for the Poisson state shown as the solid line. The 
dashed line is a plot of a%, ~ 1/m for comparison. 
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FIG. 6. Plot of &(j> versus number of atomic detections with 
Fn 2 /ni = 1,1/2 and 1/4 for the solid, dashed and dash-dotted 
curves respectively. 



FIG. 1. The conditional visibility curves for an initial num- 
ber state (n = 20), Poisson and thermal states plotted as 
dashed, solid and dash-dotted curves respectively. The visi- 
bility is dependent upon the ratio of the mean counting rates 
between the condensates (Tn^/rii). 



FIG. 2. Histogram of 500 numerically generated atomic de- 
tections plotted as circles. The solid curve is a least-squares 
fit of the form 1 + f3 cos{2nx + <f>). A thermal initial state is 
shown in (a) and an initial Poisson state is shown in (b). 
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